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We determine the phase diagram of a generahzed Shastry-Sutherland model, using a combination 
of dimer- and quadrumer-boson methods and numerical exact diagonalization techniques. Along 
special lines in the parameter space the model reduces to the standard Shastry-Sutherland model, 
the 1/5-th depleted square lattice and the two-dimensional plaquette square lattice model. We study 
the evolution of the ordered phases found in the latter two unfrustrated models under the effect of 
frustration. Furthermore we present new exact diagonalization results for the Shastry-Sutherland 
model on clusters with up to 32 sites, supporting the existence of an intermediate gapped valence 
bond crystal phase with plaquette long-ranged order. 



There has been considerable interest in recent years 
in the study of low-dimensional quantum spin systems, 
both experimentally and theoretically. Special attention 
has been devoted to two-dimensional antiferromagnetic 
systems, where quantum fluctuations and frustration al- 
low for various exotic quantum phases to compete with 
quasi-classical long-range magnetic order. 

Already in the unfrustrated regime such spin liquid 
phases occur in certain regions of the parameter space, a 
prominent example being the spin-1/2 Heisenberg model 
on the 1/5-th depleted square lattice. In this model 
two distinct spin liquid phases, well described by res- 
onating valence bond (RVB) -like states, are found along 
with an intermediate long-ranged antiferromagnetically 
ordered phase. The location of the quantum critical 
points separating these phases are known to rather high 
precision. Q 

When considering frustrated systems, ground state 
properties are less well established. A prominent example 
of recent interest, due to its relevance for the spin gap sys- 
tem SrCu2 (603)2, is the Shastry-Sutherland model 
(SSM) . 01 In this model the nearest-neighbor square lat- 
tice antiferromagnet is frustrated by additional diagonal 
interactions, arranged in a staggered pattern on alter- 
nate squares. This model retains long-range Neel-order 
for small diagonal coupling. Furthermore, it becomes an 
exact dimer valence bond solid (VBS) with singlets form- 
ing on the diagonal bonds for small axial coupling. Con- 
cerning the existence and nature of intermediate phases, 
despite numerous investigations, [8-12], a definite picture 
has not yet emerged. 



I. THE MODEL 



of the SSM. Therefore, we introduce the quadrumerized 
SSM, defined on the square lattice by the following spin- 
1/2 Heisenberg Hamiltonian 

if = S, • S, + J ^ S, • Sj +J'Y. s» • Sj. (1) 

{i.3)u {«J>n' 

Here K and J are the two inequivalent nearest neighbour 
exchange couplings, whereas J' denotes the next nearest 
neighbour alternating dimer coupling. The various cou- 
plings are displayed in Fig. 1. 
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FIG. 1. The quadrumerized Shastry-Sutherland lattice. 
A spin-1/2 degree of freedom is located on each vertex. The 
various couplings, represented by different line styles, are de- 
noted K (solid), J (dashed), and J' (dotted). The lattice is 
self-dual under the exchange (J <-> K). Arrows illustrate the 
ordered phase found for J = in the region J' k. K. 



In this paper, we present exact diagonalization stud- 
ies on the SSM that indicate the occurrence of a valence 
bond crystal (VBC) in the intermediate regime, with pla- 
quette long-range order. Furthermore, we are able to link 
this phase to a consistent phase diagram of an extension 



Note, that the lattice is self-dual under the exchange 
(J ^ K), hence only the case J < K will be considered. 
Furthermore, the standard Shastry-Sutherland model is 
recovered along J = K, and has a larger space group 
symmetry. Other lines of enhanced symmetry correspond 
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to the 1/5-th depleted square lattice (J = 0), and the 
plaquette square lattice (J' = 0), respectively. 

The paper is organized as follows: In the next sec- 
tion we present the phase diagram of the quadrumerized 
SSM as obtained from boson operator mean-field theory. 
Then in the third section we use exact numerical diago- 
nalization to study the effects of frustration in the model 
and relate the numerical results to the mean-field phase 
diagram. We concentrate on the standard SSM in the 
fourth chapter and provide evidence for a VBC interme- 
diate phase. A summary and conclusions are given in the 
final section. 

II. BOSON OPERATOR APPROACH 

We first review the numerical results obtained along 
the unfrustrated lines in the parameter space of the 
Hamiltonian (1). For the 1/5-th depleted square lat- 
tice (J = 0) there exists a plaquette RVB-like phase 
(PRVB) at small J'/K < 0.94, and a dimer RVB-like 
phase (DRVB) at large J'/K > 1.67, with an intermedi- 
ate long-range ordered antiferromagnetic phase. Q The 
classical configuration corresponding to this order is de- 
picted by the arrows in Fig. 1. 

The plaquette square lattice (J' = 0) is unfrustrated as 
well. Using stochastic series expansion quantum Monte 
Carlo simulations a quantum critical point is found at 
{J/K)c ~ 0.55, separating a PRVB spin liquid for J/K < 
{J/K)c from a gapless Neel-ordered phase. |jl^ This value 
agrees well with results from perturbation expansions. 



While quantum Monte Carlo proves powerful for 
studying the unfrustrated limits of the Hamiltonian (1), 
due to the sign problem other methods are needed, once 
frustration is present. In order to study possible insta- 
bilities of the spin liquid phases, we use standard boson 
operator mean-field theory, which is known to work on a 
qualitative level even for large frustration. Our analyt- 
ical results will also be substantiated by the numerical 
approach of the following sections. 

A. Dimer-boson approach 

Consider first the DRVB regime J' ^ J,K where the 
dimer-boson technique can be applied, [p^ In this rep- 
resentation, the spin- 1/2 degrees of freedom on each J'- 
dimer are expressed by bosonic bond operators, 

|s)=st|o)^|0,0), 

|<+)-<;|o> = |i,-i), 

|to)=4|0) = |l,0), (2) 
\t-)=tU0) = \l,+l), 



where the |S', S'^) denote the states on a given dimer. 
From the action of the spin operators Si, i = 1, 2 (denot- 
ing the two sites of a dimer) on these states the repre- 
sentation of the spin operators can be deduced 

St = ^ (4*0 + tlt^) ± ^ {tis - . (3) 

The spin commutation relations for Si, i = 1,2 are ob- 
tained when the bond operators obey bosonic statistics. 
Furthermore, the number of physical states available on 
each dimer specifies a hard-core constraint for the bosons 
on each dimer /i, 

sis^ + <L+*P.+ + *Lo*P,o + *L-ip.- = 1- (4) 

Taking the interdimcr couplings J and K into account 
and using Eq. , the Hamiltonian (1) is mapped onto 
an equivalent bosonic Hamiltonian, — + , 
containing quadratic, diagonal terms 




and quartic terms, Hj' , describing the interdimer scat- 
tering. The square lattice of dimers, with two dimers per 
unit cell, is found to reduce to a square lattice with a sin- 
gle site per unit cell in the bosonic representation. Here, 
we first imagine rotating half of the J'-dimers clockwise, 
so that all dimers align along the (l,-l)-direction (cf. Fig. 
1) . Then we take the centers of the dimers as the sites and 
use a coordinate system where ccd is along the original 
(l,-l)-direction and yo along (1,1). To proceed, we need 
to implement the constraint (|^) by means of a Holstein- 
Primakoff representation, 

= s = ^\ -tl_t+ -4*0 -tU^, 

and then decouple the quartic via a linear approxi- 
mation, similar to linear spin wave theory. The resulting 
total quadratic Hamiltonian, , is then diagonalized in 
momentum space using a generalized Bogoliubov trans- 
formation. jTrj l This approach is expected to work well 
inside the DRVB phase. A threefold degenerate spec- 
trum of triplet excitations is obtained, consistent with 
unbroken SU(2) symmetry. 



uj{kD) = ^ J' [J' + {K - J) {cos - COS ky)]. (5) 

Here, the wave vector k/j = ikx,ky)D is defined with 
respect to the dimer coordinate system (xd^uo). The 
phase boundaries of the DRVB phase are obtained from 
the instabilities of the triplet excitation spectrum, i.e. by 
a vanishing spin gap at k/j — (7r,0)£), signaling the con- 
densation of the corresponding bosons at this wave vec- 
tor, p^] Mapping back onto the original square lattice. 
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the corresponding magnetic order is obtained, character- 
izing the phase beyond the instabihty Une 2{K — J) > J' . 
In fact, this magnetic order corresponds to the long-range 
order found for the 1/5-th depleted square lattice, c.f. 
Fig. 1. 



B. Quadrumer-boson approach 

When the parameters in the Hamiltonian of Eq. (1) 
are close to another limiting case, J, J' ^ if, a sim- 
ilar approach, the quadrumer-boson technique, can be 
applied. [T^ The Hamiltonian of a single quadrumer, 
H = K[Si • S2 + S2 • S3 -I- S3 • S4 + S4 • Si], can be ex- 
pressed in terms of the total spin, S = Si + S2 + S3 + S4, 
and the total subspin on each diagonal, Sa = Si + 
S3, and Sb = S2 -f S4. The spectrum is given by 
E{\S, S'',Sa, Sb)) = K/2{S^ -S\- 3%). The lowest ly- 
ing triplet 1, = 0,±1} has a gap, A = K,to 
the ground state |0, 0, 1, 1). Since there is a further gap, 
A' = K, to the higher excitations, we attempt to obtain 
the instabilities of the PRVB phase by using a restricted 
quadrumer-boson method, omitting all the higher exci- 
tations on the quadrumers. Hence, the spin-1/2 degrees 
of freedom on each quadrumer are expressed by bosonic 
operators on the restricted Hilbert space, 



|s) -st|0> = 10,0,1,1), 
\t+)^tl\0) - |1, -1,1,1), 
|io)=4|0) = 11,0,1,1), 



(6) 



\t- 



t 



|1, +1,1,1)- 



From the action of the spin operators Si, i — 1, ...,4 in 
the restricted Hilbert space, the following representation 
can be deduced fl9fl 



< I - tit^ - 



2V2 



(4 



±^0 + to^T ) ^ 



V6 

i-iy 

V3 



t\s - sHj 



(7) 



In the restricted Hilbert space, the hard-core constraint 
is now obeyed on each quadrumer. Expressing the 
Hamiltonian (1) in terms of the quadrumer-boson opera- 
tors, a bosonic Hamiltonian, — Hq +Hf , is obtained 
with a noninteracting diagonal part 



and a quartic scattering part Hf . Here the sum extends 
over the square lattice of quadrumers formed by the K- 
bonds in Fig. 1. Following the decoupling procedure 
already used in the dimer-boson approach, the following 
threefold degenerate triplet excitation spectrum is ob- 
tained in the PRVB regime. 



w(k_f 



K — -{2 J — J') {coskx + cos ky) 



(8) 



The minimum of this spectrum is located at kp = (tt, tt) p 
when J' > 2 J, and the gap vanishes for J' > 3/A{K + 
2J), corresponding again to the order depicted in Fig. 
1. Furthermore, for J' < 2J the minimum is located 
at kp = (0,0)p, and the gap is again closed for J' < 
2J — 3/4:K. In this regime the model becomes long-range 
Neel-ordered. 

Upon comparing the ground state energies from the 
dimer- and the quadrumer-boson approach inside the 
common range of stability, we can obtain the direct first 
order transition line between the DRVB and PRVB spin 
liquid phases. |lq] 
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FIG. 2. Phase Diagram of the quadrumerized Shastry- 
Sutherland model. Ordered phases are characterized by the 
ordering wave vectors in the boson representations. Solid lines 
indicate second order transition lines, and the dashed line the 
first order transition. Also shown are points on the phase 
boundaries from quantum Monte Carlo (circles), [4,13] exact 
diagonalization (squares, and triangles), and series expansion 
(diamond). [12] Dashed-dotted lines are guides to the eye. 
The increment of the parameter scan in exact diagonalization 
was A{J/K) = 0.1 [A{J'/K) = 0.2] for the (71,71) p [(0,0)p] 
phase boundary. 

The overall phase diagram is shown in Fig. 2. The 
spin liquid phases are characterized by the correspond- 
ing RVB-like state, while we label the long-range ordered 
phases by the ordering wave vectors in the boson opera- 
tor approaches. Furthermore, solid lines represent second 
order phase transitions, whereas the dashed line indicates 
the first order transition line. 

The Neel-ordered phase (0, 0)p extends up to rather 
large frustration, with a largest extend of (J'/"^)max ~ 
1.2 along the Shastry-Sutherland line (J = K). On the 
other hand, the largest extent of the (tt, tt) p phase, for 
J'/K w 1.3 is bound by {J/K)^sx < 0.55 from exact 
diagonalization. This difference can be traced back to 
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the ratio of the number of frustrating couplings to the 
number the initial couplings, which is 1 : 4 when starting 
at J' = 0, but 2 : 3 upon starting at J = 0. 

Furthermore, from the phase diagram in Fig. 2, we 
find that the DRVB phase of the 1/5-th depleted square 
lattice is adiabatically connected to the exact dimer 
VBS phase of the standard Shastry-Sutherland model 
(the dimer VBS state fails to be an exact eigenstate for 
J ^ K). On the other hand, the DRVB phase is not adi- 
abatically connected to the PRVB, as expected on topo- 
logical grounds. Hence, we find a first order phase 
transition separating the two spin liquid phases beyond 
the regime of the (7r,7r)p phase. 

When turning to the case J > K, the phase diagram 
shown in Fig. 2 is obtained upon interchanging J and 
K, due to the invariance of the Hamiltonian in Eq. (1) 
under the exchange (J <-!■ K). Furthermore, we label the 
plaquette RVB-like phase for J > K by PRVB', since 
now singlets are predominantly formed on a different set 
of quadrumers than in the PRVB phase. 

III. EXACT DIAGONALIZATION STUDIES 

We include in Fig. 2 the positions of quantum critical 
points along the unfrustrated lines, obtained by quantum 
Monte Carlo. These compare rather well with the 

above mean-field theory. To extend the numerical anal- 
ysis into the frustrated regime we have performed exact 
diagonalization studies on clusters with iV = 8, 16, 32 
spins, using periodic boundary conditions, along various 
lines in the phase diagram. We determine the finite size 
values of the order parameter, M', defined by ||2^ 

Here et takes on the values ±1 at site i, according to the 
pattern in Fig. 1 for the (tt, tt) p phase or the standard 
Neel-order for the (0,0)p phase, respectively. Using the 
finite-size data we determine M' from the scaling law ||2^] 

+ + (10) 

For example, in Fig. 3 we show results obtained along 
the line J' = 1.3, where quantum Monte Carlo simula- 
tions at J = give a maximum moment of M' « 0.23. 

Within exact diagonalization we can reproduce this 
value, and furthermore observe a smooth decrease in M' 
upon increasing the frustration, up to a critical point at 
J/K « 0.55, where we enter into the spin liquid regime. 
Proceeding in a similar fashion we obtain the critical 
points depicted by squares in Fig. 2. Moreover, along 
the line J' = 2 J no finite order parameter was obtained 
after finite size scaling. 
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FIG. 3. Evolution of the order parameter of the (7r,7r)p 
phase upon increasing J/K for J' / K — 1.3. The inset showns 
the finite size data obtained for A'' = 8,16,32 at J/K = 
0, 0.1, 0.6, (top to bottom). 

The two spin liquid phases cannot be separated 
within exact diagonalization using space group symme- 
try. Namely, upon increasing J' /K for a constant J/K, 
the representation class of the ground state does not 
change. Nevertheless, from the approximate slopes of the 
ground state energy vs. J' /K at constant J/K in the re- 
gions of small and large J' /K respectively, we estimated 
the first order transition points indicated by triangles in 
Fig 2. 

Comparing these numerical results with the mean-field 
calculations in the last section, we conclude that the bo- 
son operator approach gives a good qualitative account of 
the phase diagram of the quadrumerized SSM. Namely, 
the characterization of the various phases and the loca- 
tion of the phase transition lines agree well with numer- 
ical results. 

IV. SHASTRY-SUTHERLAND MODEL 

In the quadrumer-boson approach we find a finite win- 
dow on the Shastry-Sutherland line (J — K), where both 
plaquette spin liquid phases, PRVB and PRVB', come ar- 
bitrary close to the J = K line. This already indicates 
an intermediate phase in the standard SSM between the 
Neel-ordered phase and the dimer VBS phase. Similar 
conclusions were obtained in a field theoretical study of 
a generalized SSM, which does not break the symmetry 
needed by the dimer VBS state to be an exact eigen- 
state. jll|] We now focus on this intermediate regime in 
the standard SSM. 

There has recently been considerable interest in the 
nature of the intermediate phase. In the large spin, clas- 
sical limit the system retains Neel-order for J' / J < 1 and 
is helically ordered otherwise, with a twist between next 
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nearest neighbour spins oi q = arccos(— J/ J'). Using 
Schwinger boson mean-field theory, Albrecht and Mila 
predicted the existence of a hehcal phase separating the 
dimer VBS and the ordered phase also for the spin- 1/2 
case, in a range 1.1 < J' / J < 1.65. ||] Field-theoretical 
studies by Chung, Marston, and Sachdcv for a general- 
ized spin-S model with Sp(2N) symmetry suggest the he- 
lical order to occur in a larger range, 1.02 < J' / J < 2.7, 
at = 1/2. [|lO| Furthermore, this approach predicts a 
phase with plaquette order in the extreme quantum limit, 
1/S > 5. Indeed, using series expansions around the 
plaquette limit of Eq. (1), Koga and Kawakami found 
the intermediate phase in a range 1.16 < J' /J < 1.48 
to be adiabatically connected to the PRVB phase. ||] 
However, extended series expansions by Weihong, Oit- 
maa, and Hamer lead to different conclusions. [|2| They 
suggest the PRVB phase to become unstable before the 
Shastry-Sutherland line is reached and found a columnar 
dimer phase to be a possible candidate for the interme- 
diate phase. 
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FIG. 4. Ground state energy per site of the Shas- 
try-Sutherland model from exact diagonalization of clusters 
with A'' = 16 (squares), 20 (triangles), 24 (diamonds), and 32 
(circles) sites. The solid line is a guide to the eye, and the 
dashed line is the exact dimer VBS state energy. 

Here, we perform exact diagonalization studies on clus- 
ters with up to 32 sites, significantly beyond the largest 
cluster sizes studied so far (24 sites in Ref. [Q). In Fig. 
4 we show the ground state energy of the SSM close to 
the transition into the dimer VBS state for various clus- 
ter sizes. The solid line shows our results for the 32 site 
cluster, whereas other sizes are represented by symbols. 
In addition, the energy of the dimer VBS state is shown 
by the dashed line. Due to the different symmetry of the 
various clusters, the ground state enery per site does not 
show a monotonous finite size dependence. Nevertheless, 
independent of the cluster size we find the system to be 
within the dimer VBS for J'/ J > 1.5, consistent with the 



upper bound for the intermediate phase given by Koga 
and Kawakami, but significantly below the values from 
the Schwinger boson, and Sp(2N) theory. From finite size 
analysis of the Neel-order parameter, we estimate a upper 
bound for the ordered phase of {J' /J)max < 1-4, consis- 
tent with the series expansion value (J'/ J) max — 1.2 ±0.1 
I p^ . More interestingly, for the largest cluster we find a 
characteristic change in the curvature of the ground state 
energy, well before entering the dimer VBS. Hence, we 
conclude that characteristic features of the intermediate 
phase could be retrieved from clusters with TV = 32 sites 
in a range 1.425 < J' /J < 1.475. 

In this regime we find the spin-spin, or 2-point correla- 
tion function to decrease rapidly with distance, indicating 
the absence of antiferromagnetic order. 
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TABLE I. Dimer-dimer correlations C''(l,2;fc,/) in the 
ground state of the Shastry-Sutherland model on the 32 sites 
cluster at J' / J = 1.45. The labeling of the sites is shown in 
Fig. 5. 
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In order to test against the various proposed ground 
states, we measure the dimer-dimer, or 4-point correla- 
tion functions, 

C\t,j,;k,l) = {S,-S, Sfe.SO-(S, •S,)(Sfc.S,), (11) 

on the iV = 32 lattice in the above interaction range. In 
particular, we fix = (0,16) and extend {k,l) over 

all inequivalent J bonds. The values obtained in the 
ground state for J'/J = 1.45 are displayed in Table I, 
and illustrated in Fig. 5, which also shows the labeling 
of the sites. 




FIG. 5. Dimer-dimer correlations in the ground state 
of the Shastry-Sutherland model on the 32 sites cluster at 
J'/J = 1.45. The reference bond is the bond (0,16). Posi- 
tive (negative) correlations are drawn as full (dashed) lines. 
The thickness of a line is proportional to the strength of the 
correlation. Short diagonal lines indicate the position of the 
J'-dimer bonds. 

We obtain a clear signal in the dimer-dimer correla- 
tions that extends throughout the whole cluster, with 
a finite asymptotic value approximately reached for the 
larger dimer-dimer distances on the cluster. In the spatial 
distribution we furthermore observe periodic oscillations, 
which reflect an underlying order of quadrumer-singlets, 
formed predominately on void squares (i.e. those squares 
not containing a diagonal bond). 

For the SSM there are two equivalent configurations 
with quadrumer-singlet coverings residing on the two dif- 
ferent subsets of void squares (formed by the J or K 
bonds in Fig. 1, respectively). Hence, if the plaquette- 
like order in the 4-spin correlation function survives 
quantum fiuctuations, indicative of a plaquette VBC, a 
two-fold degenerate ground state manifold will emerge in 
the thermodynamic limit. On a finite lattice this degen- 



eracy is lifted, but a low-lying singlet state well inside the 
triplet gap and only slightly above the ground state en- 
ergy is expected. We can obtain the quantum numbers of 
this singlet state from the following symmetry considera- 
tions: The Shastry-Sutherland lattice has a p4mm space 
group symmetry, Q and the ground state of the 32 site 
cluster has momentum k = (0,0), and is invariant under 
the 7r/2 rotations about the center of any void square (i.e. 
s-wave like). Furthermore, the two equivalent configura- 
tions of quadrumer-singlet coverings are related by the 
reflections about either diagonal dimer axes, but are in- 
variant under the tt/2 rotations and lattice translations. 
Namely, they both reside already inside a single unit cell 
of the Shastry-Sutherland lattice. 
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FIG. 6. Ground state energy, and low-lying singlet and 
triplet excitations in the k = (0, 0) sector of the Shas- 
try-Sutherland model on the 32 sites cluster. Open symbols 
represent singlet states, full symbols triplet states. Circles 
denote states with eigenvalues ~ 1 (s-wave), squares de- 
note R„/2 = — 1 (d-wave), and triangles R„/2 ~ H (two-fold 
degenerate) . The solid line is a guide to the eye for the ground 
state, the dotted to the lowest excited s-wave state, and the 
dashed line the lowest d-wave state. The Dashed-dotted line 
is the exact dimer VBS state energy. 

Hence, a low-lying s-wave symmetric singlet state with 
momentum (0,0) is expected to be included in the spec- 
trum of the 32 site cluster in the regime of the plaquette 
VBS phase. This state should furthermore show similar 
dimer-dimer correlations as the absolute ground state. 

In Fig. 6 we plot the ground state energy along with 
those of the lowest excited singlet and triplet states for 
the 32 site cluster in the zero momentum sector. We spec- 
ify the transformation properties of the various states 
under the 7r/2 rotation by the eigenvalue i?7r/2 = 1 (s- 
wave), —1 (d-wave), or ±i (two-fold degenerate). See 
the caption of Fig. 6 for a detailed account on the var- 
ious used symbols. In the regime where we expect ev- 
idence for an intermediate state we indeed find various 
low-lying singlet states well inside the rather large triplet 
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gap. Moreover, there are two singlet states with ener- 
gies rather close to the ground state energy, one being 
a s-wave, and the other a d-wave state with respect to 
the 7r/2 rotations. We furthermore calculated the dimer- 
dimer correlations for both states and find for the s-wave 
state a similar signal as for the absolute ground state. 
On the other hand, the d-wave state does not show any 
pattern in the dimer-dimer correlations. This state does 
not seem relevant for the ground state of the system in 
the thermodynamic limit. Presumably it is related to 
low-lying excitations in the dimer VBS phase. p^j2^ 

From the above we conclude that the low-lying s-wave 
singlet state, having the right quantum numbers and 
dimer-dimer correlation, will become degenerate with the 
ground state upon increasing the cluster size to the in- 
finite lattice. Both states then form the two-fold degen- 
erate ground state manifold of a plaquette VBC in the 
thermodynamic limit. In Fig. 7 a pictorial description of 
this ground state manifold is given. 
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FIG. 7. Pictorial illustration of the variational ground 
state manifold of the Shastry-Sutherland model in the pla- 
quette VBC phase. Thick lines indicate four spins involved in 
a quadrumer-singlet. 

Furthermore, this ground state manifold is invariant 
with respect to lattice translations and the 7r/2 rota- 
tions about the centers of the void squares, but each state 
spontaneously breaks the reflection symmetry about the 
diagonal axes along the dimer directions. 

From our numerical results we can also exclude a 
columnar dimer state, which would be four-fold degen- 
erate and furthermore show a different pattern in the 
dimer-dimer correlations than depicted in Fig. 5. ]2^ ] 

Further evidence for the relevance of the plaquette 
VBC in the intermediate regime of the SSM can be 
drawn from analogous results of recent studies on the 
spin-1/2 Heisenberg model on the checkerboard lattice. 
P^ , ^ In this model, diagonal bonds are again organized 
on a square lattice in a pattern as to leave half of the 
squares void. Namely, the underlying lattice is obtained 
by adding in an additional diagonal bond on each square 
on which a dimer bond is located in Fig. 1. Also in 
the checkerboard lattice quadrumer-singlets form on the 
void squares, resulting in a two-fold degenerate plaque- 
tte VBC ground state manifold with broken space group 



symmetry. The structure leading to frustration is rather 
similar in both models, and the system tries to minimize 
frustration by forming plaquette valence bonds on the 
void squares. Due to the homogeneous axial couplings, 
this can be accomplished only upon spontaneously break- 
ing the symmetry inherited from the underlying lattice. 



V. CONCLUSIONS 

In conclusion, we have studied the phase diagram of the 
quadrumerized Shastry-Sutherland model. Using bond- 
operator methods and exact numerical diagonalization its 
phase diagram was established, which links the various 
results available for special limiting cases of the model. 
The antiferromagnetically ordered phase of the 1 / 5-th de- 
pleted square lattice model is destroyed by modest frus- 
tration, whereas the Neel-ordered phase extends up to 
rather large frustration. There is a first order transition 
line separating the different spin liquid phases, PRVB 
and DRVB, beyond the ordered phase. Furthermore, the 
DRVB is adiabatically connected to the dimer VBS phase 
of the SSM. 

For the standard SSM there exists a finite region 
around J'/ J = 1.45 where the system becomes a plaque- 
tte VBC with spontaneously broken space group sym- 
metry, and a two-fold degenerate ground state manifold. 
Perturbing away from the Shastry-Sutherland line, i.e. 
for J 7^ if, the symmetry is broken explicitly, and the 
system favors a unique ground state, namely the PRVB 
for J < K, and the PRVB' for J > K. Furthermore, 
upon varying J' along the Shastry-Sutherland line, a first 
order transition leads to the dimer VBS and a second or- 
der transition to the Neel-ordered phase. Within our an- 
alytical and numerical studies we did not find indications 
for further phases in the (quadrumerized) SSM. 

Stabilization of a plaquette VBC phase in both the 
Shastry-Sutherland and the checkerboard lattice model 
agrees with the generic structure of the underlying frus- 
trated lattice. However, in the SSM the VBC is unsta- 
ble towards the dimer VBS upon increasing the diagonal 
coupling. For the checkerboard lattice model the range 
of the VBC phase is still unkown, and remains for further 
studies. 

We acknowledge fruitful discussions with Stephan 
Haas, Andreas Honecker, and Bruce Normand. This 
work has been financially supported by the Swiss Na- 
tionalfonds. 

Note added.- After completion of this work Akihisa 
Koga, and Norio Kawakami pointed us to Ref . , where 
the phase diagram of the Hamiltonian in Eq. (1) was 
studied using series expansion methods. Their results are 
in perfect agreement with our calculations in sections II 
and HI. However in the VBC nature of the inter- 
mediate phase in the Shastry-Sutherland model was not 
noticed. 
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